A study of neutron-deuteron scattering in configuration space 
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A new computational method for solving the configuration-space Faddeev equations for 
the breakup scattering problem [J has been applied to nd scattering both below and above 
the two-body threshold. To perform numerical calculations for a general nuclear potential 
and with arbitrary number of partial waves retained, we use the approach proposed in [2] . 
Calculations of the inelasticities and phase-shifts, as well as cross-sections and analyzing 
powers for elastic nd scattering at Ei a b =14.1 MeV were made with the charge independent 
AV14 potential. The results are compared with those of other authors and experimental 
data. 



1. Formalism 

The Faddeev equations projected onto the MGL basis |2] are : 
E + -{dl + d 2 y ) - v%(x,y)]<f>*» s °> M °(x,y) 
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<Z>^ so > Mo (x,y) + I' duJ2g 0J (x,y,u)<Z>^ M °^' ^ 
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Here Greek subindexes denote state quantum numbers: a = {l,a, J,s,X,t}, where I, a, 
J and t are the orbital, spin, total angular momenta and isospin of a pair of nucleons, A 
is the orbital momentum of the third nucleon relative to the cm. of a pair nucleons, and 
s is the total "spin"(s = 1/2 + J). M = A + sis the three-particle angular momentum. 
The geometrical function gp 7 (x,y,u) is the representative of the permutation operator 
P + + P in MGL basis (details see in [2]). In Eq. (1) v£ is the centrifugal potential, and 
nucleon-nucleon potentials ) = SxySss'Saa'Sjj'Vf/ , where v^, are the potential 

representatives in the two-body basis 3^ Jz (x). 
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The asymptotic conditions for the radial parts are 
^ s ° M °(x,y) = 8 XXo 8 SSQ 8 (Tl 8j 1 Uqy)H^) + i<° XoSQ h{ +) (qy)Mx) + O^ 1 ) (2) 
and 

^ M (x,y) = A^ oSQ e - 7T + 0(X-^), X 2 = x 2 + y 2 . (3) 

In these formulae j\(qy) is the regularized spherical Bessel function, ipi(x) is the deuteron 
wave function (1=0,2), h x + \qy) is the regularized Hankel function, a^° x and A^ XoSQ are 
the elastic and breakup amplitudes, respectively, q an E are the cm momentum and final 
energy. 

Table 1 



n — d elastic amplitudes. 



a \s,X Q SQ 


Ref. |3] 


present 


Ref. 


present 




Elab=l MeV 




Elab=3 MeV 




1/2+ 
a i/20,l/20 

1/2+ 
a 3/22,3/22 


-0.291 + i 0.093 
-0.0175+ i 0.0003 


-0.2903+ i 0.093 
-0.0173+ i 0.0003 


-0.4691+ i 0.3272 
-0.0682+ i 0.0048 


-0.4679+ i 0.324 
-0.0616+ i 0.004 


1/2" 
a i/21,l/21 

1/2" 
a 3/21,3/21 


-0.0718+ i 0.0055 
0.2085+ i 0.0459 


-0.0699+ i 0.0052 
0.208 + i 0.0457 


-0.122 + i 0.0198 
0.3749+ i 0.176 


-0.120 + i 0.0191 
0.3739+ i 0.1743 


3/2+ 
a 3/20,3/20 

3/2+ 
a i/22,l/22 

3/2+ 
a 3/22,3/22 


-0.4985+ i 0.5383 
0.0101+ i 0.0001 
-0.0187+ i 0.0004 


-0.4985+ i 0.5382 
0.0108- i 0.0001 
-0.0181+ i 0.0004 


-0.3145+ i 0.888 
0.0420+ i 0.0018 
-0.0735+ i 0.006 


-0.3153+ i 0.8874 
0.0437- i 0.0019 
-0.0721+ i 0.0057 


3/2- 
a i/21,l/21 

3/2- 
a 3/21,3/21 

3/2- 
a 3/23,3/23 


-0.0718+ i 0.0052 
0.2392+ i 0.061 
0.0021+ i 7.1E-6 


-0.0716+ i 0.0052 
0.2401- i 0.0615 
0.0023+ i 7.6E-6 


-0.1227+ i 0.016 
0.396 + i 0.1959 
0.0158+ i 0.0003 


-0.1212+ i 0.0157 
0.3974- i 0.1978 
0.017 +i 0.0004 



2. Results 

In the present run we used the charge independent AV14 potential and included all 
angular momenta of subsystems, I and A < 4, total angular momentum of a pair nucleons 
J < 3, and total three-body angular momentum Mo up to 9/2. The cutoff radius was 
taken 80 fm. In Table^our results on the nd elastic amplitudes are given. For comparison, 
accurate results of Bochum and Pisa group are also presented. 

In Figure ^ the elastic differential cross-section is shown together with that of Bochum 
group [3] and the experimental data of |4|5j . Agreement is good except for angles in the 
vicinity of 180 degrees. In Figures 121 and El vector analyzing powers are shown along with 
the experimental data. 
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Figure 1. nd elastic differential cross section at Ei a b=14.1 MeV. The solid line is our 
results. The dashed one corresponds to the Bochum group prediction [3]. The data are 
from |lj (open circles) and [Sj (open diamonds). 




Figure 2. The vector analyzing power Ay at 

Ei ah = 14.1 MeV. The solid line is our results. Figure 3. The deuteron analyzing power iT n 
The dashed one is the Bonn B prediction with at Ei ab = 14.1 MeV. The solid line is our results. 
J <3 of 6|. The data are from jjj. The pd data at Ei a b=15.0 MeV are from [Hj. 

As one can see agreement with experimental data is enough good in both cases. How- 
ever one to take into account our limitations on the maximal values for M, I and A. 
Convergence with these values is presently under study. 

To check accuracy of our calculations we applied the optical theorem. Presenting it in 
the form Ima = R.H.S, where the expression for R.H.S may be found in we get the 
results shown in Table El 
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Table 2 

The optical theorem results for E lab = 14.1 MeV and M£ = l/2+,3/2 



s 0> ^0 


t < 2+ 

AoSo,A()S() 


R.H.S. 


s 0; \ 


Imaf° c x „ 


R.H.S. 


1/2,0 


0.7172 


0.7182 


3/2, 


0.8942 


0.8971 


3/2,2 


0.0258 


0.0263 


1/2,2 


0.0516 


0.0522 








3/2,2 


0.0359 


0.0363 



Presenting the total cross-section a to t as a sum of its elastic a e i and inelastic a in parts 
we obtained values for them shown in Table El 

Table 3 

Cross-sections of nd breakup scattering for Ei a b =14.1 MeV in [mb] 



M 


Otot 


0"eZ 




1/2 


183.25 


132.47 


50.78 


3/2 


584.82 


471.14 


113.68 


5/2 


786.61 


637.14 


149.47 


7/2 


804.25 


648.20 


156.05 




824. ± 10 


694. ± 42 


130. ± 43 




807. |] 




172. ± 12 HO] 



3. Conclusion 

1. The proposed method of the solution of the Faddeev equations in configuration space 
by the direct integration of the set of partial differential equations is fully feasible and 
reliable, and give results of the same precision as obtained by different methods. 

2. This shows that the method can be generalized to include the Coulomb interaction in 
the pd scattering below and above the deuteron threshold, where others methods experi- 
ence considerable difficulties especially for the breakup channel. 
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